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Synthetic and biological living polymers are self-assembling chains whose chain length distribu- 
tions (CLDs) are dynamic. We show these dynamics are ultrasensitive: even a small perturbation 
(e.g. temperature jump) non- linearly distorts the CLD, eliminating or massively augmenting short 
chains. The origin is fast relaxation of mass variables (mean chain length, monomer concentration) 
which perturbs CLD shape variables before these can relax via slow chain growth rate fluctua- 
tions. Viscosity relaxation predictions agree with experiments on the best-studied synthetic system, 
Q-methylstyrene . 
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The term "Living Polymers" labels a remarkable 
and diverse family of self-assembling systems in which 
molecules or other microscopic units spontaneously ag- 
gregate into long chains by continuously adding to chain 
ends. Their technological and biological importance and 
their unusual properties as examples of soft condensed 
matter have driven a large body of experimental and the- 
oretical research^ H H 0, H H H B [3 [H [13 The 
classic synthetic example is ionic living polymerization^, 
0, 0, 0, Hi IE where monomers assemble into flexible 
polymer chains by adding to charged chain ends which re- 
main "alive" even after monomer is consumed, a property 
widely exploited commercially to synthesize high perfor- 
mance block copolymers and other novel materials'^. A 
recent variant on this theme with immense technological 
potential is living free radical polymerization(8| where 
ingenious capping-decapping schemes prevent termina- 
tion and yet permit chain propagation. In the biological 
world, actin and microtubule filaments, rapidly assem- 
bled from the proteins actin and tubulin, are essential to 
the motility and structural integrity of living cells . 

As polymeric or filamentous materials, the novel and 
distinguishing feature of these systems is that unlike in- 
ert polymers the chains are dynamic or "living" objects 
whose lengths constantly fluctuate. The polymerization 
processes are alive: when external conditions change, 
chain length distributions can respond dynamically and 
relax to a new equilibrium. This adaptability is the cru- 
cial property exploited in both synthetic and natural ap- 
plications. For example, it enables living cells to rapidly 
initiate motion or shape changes by altering cellular con- 
ditions in response to extracellular signals. 

In this letter we study this dynamical responsiveness 
theoretically. Whilst a rather clear picture of equilibrium 
properties has been established, living polymer dy- 
namics are far less well understood 0, la, S • We focus 
on the synthetic system which has received most exper- 
imental attention0 0, the anionic living polymer 
poly a-methylstyrene (PAMS) and we compare our pre- 
dictions to PAMS viscosity relaxation measurements at 
the end of this letter. Aside from its intrinsic impor- 



tance, PAMS is a model system for more complex cases 
such as biological living polymers and most of our re- 
sults are completely general. We use the term "Living 
Polymers" in the traditional sense to denote sys- 
tems where (i) the concentration of living chains, Pchains, 
is fixed for all time by the number of initiators, and (ii) 
chains grow at their ends only. Related classes include 
systems where surfactants aggregate into spherical 
or elongated ( "wormlike" ) micelles 1 121 whose number is 
not fixed, and whose dynamics 1 1 li Il3 | are very different. 

The principal conclusion of this work is that living 
polymers are ultrasensitive, i. e. highly dynamically sus- 
ceptible to small perturbations. Thus a small change in 
external conditions, inducing a small change in the equi- 
librium state well-described by a linear susceptibility, has 
nonetheless a large dynamical effect: intermediate states 
deviate strongly from equilibrium in that some observ- 
ables are perturbed in a highly non-linear manner. To 
quantify this, consider PAMS whose equilibrium chain 
length distribution (CLD) is close to the broad ex;ponen- 
tial 00 predicted by theory, (j)eqiN) ^ e^^/^'" (see 
fig. with typical mean number of monomer units 
per chain Neq in the range of a few 100 to several 
1000. A standard experimental procedure is the "T- 
jump" where temperature is suddenly changed by an 
amount AT. Consider a PAMS system with N^q = 1000 
subjected to a small decrease AT = — 5"C as measured 
by the small parameter e = — AT/Tq « 0.1, where 
To EE dT/dh\meq « 50°C for typical PAMS studies and 
TUeq is the equilibrium monomer concentration. The the- 
ory presented here predicts the onset after '-^ 3 hours of 
a drastic depletion of short chains. By 10 hours a hole 
has appeared in the CLD in the region < AT < 400 (see 
fig. ^1 and lengths N ~ 100 have virtually disappeared. 
This is despite the fact that the change in the equilib- 
rium CLD (t)eq{N), which is recovered after ^ 10'^ hours, 
is destined to be small, 5(j)eq{N) / (f)eq{N) = 0{e) for all 
N . This hole completely invalidates perturbation theory: 
there are no small perturbations (other than those so tiny 
as to be beyond typical experimental resolution) . This is 
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FIG. 1: Numerical results for response of CLD to a small 
T-jump perturbation of amplitude e = 0.1 perturbing T- 
dependent parameters v~ and in eq. Q. Final equi- 
librium exponential (curve f, Neq — 1366, 8 = 0.37) is close 
to initial equilibrium exponential (curve i). Along its path 
i f the CLD deviates non-linearly from both. Times 
in units of v~ (for PAMS Q, V ^ 0.1s). Timescales: 
f ^ 400,rfast « 3640, rmi x 1.32 x 10^T,low ~ 1.86 x 10^ 
Curve a, t = 3640: a hole has appeared at small A''. Curve b, 
t = 3 X 10": hole-filling. Curve c, t = 2.3 x 10^: final shape 
relaxation towards equilibrium, curve f. Inset: Relaxation of 
fast and slow variables. Arrows from left to right indicate 
t* ,r(s.st,Tmi,Taiow Fast variables Smt/m^q = —9 5Nt/Neq 
(solid) decay exponentially after rfast, and are then enslaved 
to the slow variable —54>t{0) (dashed). Note initial nonlinear 
increase of —S(j)t{Q) after t* . Theory predicts 54>t{Q) ~ t'^^"^ 
during Tm <t< rsio„ and (5<^t(0) ~ e"*/^=i™ for t > TsIov,- 

the essence of ultrasensitivity. 

We consider living polymers, such as the PAMS sys- 
tems studied by Greer and coworkers P, 0, IE ' where 
monomers spontaneously polymerize below a ceiling tem- 
perature Tp, see fig. 121 The dynamics of (j>t.Ni the 
number of chains of length N, are ipt^N = v^4't.N-i ~ 
{v'^ + v~)(j>t,N + v~(j)t,N+i where v'^ = k^m and 
are monomer addition and dissociation rates from chain 
endsfisj {k'^ is the addition rate constant). Rearranging 
terms and taking the continuous limit immediately gives 

dd^tldt = -vtd(l>t/dN + Dtd'^c^t/dN^ , 
t;t0t(O)- A[50t/5iV]Ar=o = O , (1) 

where a zero current boundary condition applies at iV = 
0. Note a second derivative term naturally emerges, with 
a "diffusivity" coefficient Dt = (i;+ + u^)/2, representing 
fluctuations about the average growth rate or "velocity" 
Vt = v'^ — v~ . Both vt and Dt depend on (pt (see eq. 
(fTj) and following remarks). Setting d(j)t/dt = gives 
the equilibrium exponential CLD with Neq = —Dcq/voq. 
Since Neq ^ 1, it follows that in equilibrium w+ « , 




FIG. 2: A living polymer is an initiator (square) plus A'' 
monomers. Monomers can add (subtract) from its live end 
(filled circle) at rates (^~)- The initiators of refs. 0,0 are 
bifunctional (one interior initiator plus two live ends). 

the diffusivity is D^q « w~ and the velocity has a small 
negative value relative to the characteristic scale 

Veq = -V'/Neq . (2) 

This is just sufficient to negate diffusive broadening 
which would otherwise smear out the equilibrium MWD 
of width Neq after time t^iow = Nl^/D^q, i. e. , Woq7"siow = 
—Neq- Note that since the average velocity over all chains 
must vanish in equilibrium (there can be no net dissoci- 
ation to the monomer pool) the > chains must have 
a small negative velocity to balance the unique = 
chains (consisting of an initiator only, fig. \^ which can- 
not depolymerize and so have a positive velocity. 

A clue as to the origin of ultrasensitivity is already 
apparent in this equilibrium situation. The significance 
of eq. ((2Jl is that equilibrium is an extremely delicate 
balance sustained by a tiny negative velocity easily over- 
whelmed by even a very small perturbation. Consider an 
equilibrium system suddenly subjected to a small nega- 
tive T-jump AT, after which the current monomer con- 
centration m exceeds its equilibrium value for the new 
temperature by 5ra w —{AT/To)meq = erueq {SX de- 
notes the deviation of any property X from its equilib- 
rium value, Xcq, which will eventually be attained after 
relaxation). This generates a velocity boost 6v — k^Sm. 
Since k^rueq ~ v~ to order 1/Neq, we have 

Sv « et)~ . (3) 

Comparing eqs. (0) and Q one sees that despite the 
smallness of the perturbation, e ^ 1, the delicate bal- 
ance of velocity and diffusion which characterizes equi- 
librium is destroyed since the boost Sv greatly exceeds 
the equilibrium velocity Vcq. This holds for any e above 
a tiny threshold 1/Neq. Immediately after the pertur- 
bation, diffusive broadening ~ {Deqt)^^^ is still stronger 
than coherent chain growth ~ ev~t. But for times larger 
than t* « Deq/€^v~ chains grow coherently, i. e. N will 
increase until the excess mass 6m held by the monomer 
reservoir has been transferred to the polymer system and 
the velocity boost has decayed. In the process the mean 
chain length increases by an amount —6N determined by 
conservation of total concentration of monomers wtotai: 

m +Pchains ^ = n^total , 5N / Neq ^ -f. 6^'^ .(4) 

The parameter 9 = (mtotai - meq)/meq « (Tp -T)/Tq 
measures distance into the polymerization regime, and 
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will be taken as order unity here. Thus, the entire CLD 
translates uniformly by « eN^q (fig. ^ in time 

Ti^st^SN/{€v-)=N,J{ev-) . (5) 

This leaves behind a hole in the CLD: chains with lengths 
less than eN^q have disappeared. 

For times beyond Tfast , mass transfer is essentially com- 
plete and TO, V and N are very nearly relaxed. These are 
the fast variables. The process of CLD shape relaxation 
needs much more time. It relies on the far slower diffusive 
process of incoherent reshuffling of monomers between 
chains. The time Tgu needed to fill the hole is simply the 
diffusion time corresponding to the hole width, 

Tfiii = . (6) 

The last process is global CLD shape relaxation on the 
scale Neq, requiring a diffusion time Tgiow — N^^/Dcq. 
These events are depicted in fig. ^ 

In summary, relaxation to the new equilibrium state 
involves 3 distinct episodes: (i) coherent chain growth 
for < i < Tfast during which fast mass variables re- 
lax, (ii) hole-filling, i. e. recovery of short chains, during 
Tfast < t < Tfiii and (iii) global diffusive relaxation for 
t > Tfiii during which slow shape variables relax on a 
timescale Tgiow. It is this conflict of timescales which is 
the origin of the non-linear hole produced by (i) : slow dif- 
fusive shape equilibration simply cannot keep pace with 
the rapid deformation produced by mass transfer. 

It is interesting to compare this with spherical micelle 
aggregation which involves two distinct relaxation pro- 
cesses llj. The first entails fast mass exchange between 
monomers and nearly monodisperse micelles whose num- 
ber remains essentially fixed due to large aggregate nucle- 
ation/dissociation barriers. This is similar to process (i) 
above. However, a crucial difference is that the number 
of aggregates can ultimately change and the timescale 
for this process in consequence depends inversely on to- 
tal monomer concentration (whereas Tfast is independent 
of TOtotai). During a second, much slower process, the 
number of micelles re-equilibrates. This process has no 
analogue for the living polymers we study. 

Let us now briefly outline formal calculations justify- 
ing these arguments, starting from eq. Its steady 
state solution is the equilibrium Flory-Schultz CLD, 

(^,qiN) = e-^/^^V^eg with N,q « mtotal0/[Pchains(l + 

0)]. Multiplying eq. Q by and integrating one has 

it = -Wi/Tfast - A0t(O)/Tfast , (7) 

after using eqs. Q and ©. Given the CLD 0t, this 
relationship determines the time-dependent velocity vt 
and thence diffusivity Dt — Vt/2 + v~ . The technical 
difficulty is that both vt and Dt depend non-locally on 
i. e. eq. is a non-linear and non-local system. 
Consider a negative T-jump perturbation inducing, as 
discussed, a velocity ?;o ~ ev~ at t = 0. Now, to order e, 



Dt K, v~ Ki Dcq- Thus in the velocity kinetics, eq. lO, 
the Vt term is initially much greater than the Dt term 
(provided e > 1/Neq). Hence vt ~ ^y- g-t/^s^nt jj^ ^j^g 
CLD evolution kinetics eq. Q in which, for t > t* , the 
coherent term wins and the CLD translates along the N 
axis, (t)t{N) K, (I)q[N — Vt). Its displacement converges 
for t :s> Tfast to eNeq/0. In other words the CLD trans- 
lates and then halts on a timescale Tfast, leaving a hole 
of size « eN in its wake (see fig. More precisely, the 
trailing edge broadens by (Doq^)^^^ and produces an ex- 
ponentially small amplitude at the origin, (j)t{0) ~ e~*/* . 

This concludes episode (i). The linearly related fast 
variables v, m and N have the same decay kinetics and 
for t 3> Tfast have all relaxed. Since the velocity is ex- 
ponentially small, the CLD kinetics eq. Q now describe 
essentially pure diffusion with refiecting boundary condi- 
tions. This is episode (ii): on the timescale Tgn, diffusion 
fills the hole (see fig. and replenishes the amplitude 
at the origin, 0t(O) « {t/Tfin)^^^e-'^^''''/y Neq- This re- 
covery of 0t(O) and the decay of vt imply that at time 
Tqs ~ (Tfast Tfiii)^/^ thc two RHS terms in the velocity 
dynamics eq. ^ become equal. Self-consistently, one 
finds that thereafter they remain matched, i. e. Vt is much 
smaller than either of these terms. Thus from this time 
onwards, the fast velocity variable (intrinsic timescale 
Tfast) evolves quasi- statically, enslaved to the slow vari- 
able (/'t(O) (intrinsic timescale Tgiow) according to 

vt- -v-(f>t{0) (i>Tqs) . (8) 

It follows that the velocity now undergoes a recovery 
as the hole fills up. By Tun its magnitude is of order 
the equilibrium value and can thus compete with diffu- 
sion. This heralds the onset of episode (iii): we return to 
the basic dynamics, eq. Q), with velocity and diffusion 
terms now of equal importance. For thc first time, rel- 
ative perturbations of all quantities are now small, and 
we can apply standard perturbation theory. We find the 
= chains recover as (50t(O) ^ t^^/^ up to TsIow and 
^ g-t/i-aiow foj. f > Tsiow. During these very late stages 
the fast variables, which long ago decayed close to equi- 
librium, are fine-tuned to their true equilibrium values, 
following 0t(O) quasi-statically according to eq. (|HJl. 

Note we implicitly assumed sufficient time for the hole 
to develop before mass transfer is complete, i. e. Tfast > 

~ 1 /2 

t* . This is true provided e > Cc = l/Neq . For very 
small perturbations e < Ec the relative depth of the hole 
though much greater than e is less than unity. Finally, we 
have studied AT < 0. The response to a positive T-jump 
is similar, but instead of a hole a large peak ^ Neqe-"^ 
develops for small A^, decaying after Tfin. 

We tested our theory by numerical integration of cou- 
pled ODEs describing the monomer-polymer dynamics. 
Fig. n shows the predicted deep hole developing after a 
small negative T-jump. Long time enslavement of fast 
variables to slow ones is clearly demonstrated (inset). 
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FIG. 3: Viscosity relaxation measurements by Ruiz-Garcia 
and Greer 0] following small T-jumps (|AT| « 1°K) on 
PAMS in tetrahydrofuran initiated by sodium naphthalide. 
Reproduced from fig. 1. of ref. 0. Circles (triangles) 
refer to positive (negative) AT. mtotai = 0.29 gm/cm'^, 
Pchains ~ 2.7 X 10~* gm/cm'^ (semi-dilute conditions). 

We conclude by discussing our results and comparing 
to experiment. We found that, following a perturba- 
tion, living polymers relax by adjusting (a) their total 
mass, proportional to first moment N and (b) the shape 
of their CLD. Mass variables are fast (relaxation time 
Tfast ~ N) because monomer-polymer mass transfer is a 
coherent process. Coherent chain growth, however, can- 
not affect CLD shape whose relaxation therefore relies 
on slow diffusion-like fluctuations in chain growth rates 
(relaxation time Tgiow ~ iV^.) These timescales are typi- 
cally separated by 2 or 3 orders of magnitude. In fig. O 
we reproduce viscosity relaxation measurements by Ruiz- 
Garcia and Greer 7] after small positive and negative 
temperature jump perturbations of the living polymer 
PAMS at 6 different temperatures. Now generally we ex- 
pect viscosity rj ~ C(j,N'^ where 7 is a characteristic expo- 
nent (31 and C0 depends on the entire CLD, 4>{N), includ- 
ing shape properties. Our theory thus predicts a fast ini- 
tial relaxation of 77 in a time Tfast followed by a very slow 
relaxation in Tgiow Now Pchains ~ 2.7 x 10"'* gm/cm'^ 
while fc+ « 0.2M~* sec~* was measured to vary by 
« 10% over this temperature range. Thus, rewriting 
Tfast = 1 /pchains^"*" , we predict a constant relaxation time 
(±10%) for all 12 measurements, Tfast ~ 2000 seconds. 
This agrees very closely with experiment (see fig. de- 
spite the fact that the observed viscosities varied by an 
order of magnitude and the mean chain lengths are esti- 
mated to vary from almost zero up to several thousand 
units at the lowest temperature. Note Tgiow = ■'"f^st''^^^ ~ 
1 week-1 month, after estimating « 0.5 and using 
v~ =0.1 sec~^. The second shape-derived relaxation of 
77 is therefore unobservably long for this experiment. 



Very large temperature quenches were studied in ref. 
0. TV relaxed faster than the second moment, though it 
was unclear if m and N were coupled as required by mass 
conservation, possibly due to ionic aggregation IO3. The 
authors of ref. studied large perturbations theoretically, 
identifying the timescales Tfast and Tgiow ■ Their analytical 
results were valid for short times and very small 9. 

In this letter we showed that even a small perturbation 
leads to a non-linear dynamical response. This ultrasen- 
sitivity is due to the inability of slow CLD shape variables 
to keep pace with the fast relaxation of N which entails 
simple translation of the CLD leaving a hole or peak at 
small N. Of all slow variables, the most sensitive is 4>t{0), 
the number of iV = chains, which becomes exponen- 
tially small or massively enhanced. These free initiators 
also govern the very late fine-tuning of the fast variables 
N and to, the free monomer concentration. Physically, 
this is because monomers add to chain ends only. But 
all ends are identical except for the = chains which 
cannot shed monomers. Thus 0t(O) is the only dynamic 
polymer property featuring in the kinetics of m. Given 
its central role, we propose measurement of the number 
of free initiators by spectroscopic or other methods as a 
revealing probe of ultrasensitivity. 

This work was supported by the Petroleum Research 
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